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Abstract 

We consider a stochastic partial differential equation (SPDE) on a lattice 

d t X = (A - m 2 )X - XX P + n 

where r] is a space-time Levy noise. A perturbative (in the sense of formal power series) strong 
solution is given by a tree expansion, whereas the correlation functions of the solution are given 
by a perturbative expansion with coefficients that are represented as sums over a certain class 
of graphs, called Parisi-Wu graphs. The perturbative expansion of the truncated (connected) 
correlation functions is obtained via a Linked Cluster Theorem as a sums over connected graphs 
only. The moments of the stationary solution can be calculated as well. In all these solutions 
the cumulants of the single site distribution of the noise enter as multiplicative constants. To 
determine them, e.g. by comparison with a empirical correlation function, one can fit these 
constants (e.g. by the methods of least squares) and thereby one (approximately) determines 
law of the noise. 
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1 Introduction and overview 

In this work we give a sufficiently general perturbative solution to the Stochastic Partial Differ- 
ential Equation (SPDE) 

= AX(t, x) - m 2 X{t, x) - XXP(t, x) + jy(t, x) , A > 0, m > 0,p G N (1) 



X(0, x) = /(x), (t, x) e]0, oo[xL 5 

where rj(t,x) is a general space-time noise of Levy type with moments of arbitrary order such 
that such that one can extract information on the distribution of ij by comparison with empirical 
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data. To avoid so-called ultra-violet (UV) divergences, we consider this equation on a spacial 
lattice = {5z , z € Z d } with lattice spacing 5 > 0. 

Nonlinear SPDEs driven by non Gaussian noise have been discussed both from the physical 
[141 Section 4.2 and Refs.] and the mathematical side, see e.g. 013]. But so far there seems to be 
no recipe to find out which kind of noise - given the general structure of the equation, including, 
say, the numerical values of all coefficients the right hand side - explains a set of empirical data 
best. It is the aim of this article to give such a recipe by calculating the correlation (moment) 
functions of the solution X perturbatively. To do so, we generalize the classical Feynman graph 
approach of Parisi and Wu ^S] ° n the Gaussian stochastic quantization equation to the general 
Levy case, using a generalized class of graphs, henceforth called Parisi- Wu graphs. Statistically 
relevant quantities as expectation, variance and higher order cumulants of the noise enter in 
the expansion as simple constants - like additional coupling constants - and can therefore be 
extracted rather easily by comparison with the data. 

The restriction to equation is not essential at all for the combinatorial graph calculus 
that we develop. The following generalizations are rather obvious: i) the nonlinearity X p can be 
replaced by a polynomial, e.g. the Mexican hat potential 2X(X 2 — a) in equations of reaction- 
diffusion type; ii) nonlinearities include derivatives, like |VXj 2 in the KPZ-equation hi) 
in addition the linear operator content on the right hand side is of higher order and the noise is 
of derivative type as in the non-linear description of nuclear beam epitaxy, where the right hand 
side of (0) is -KA 2 X - AA|VX| 2 + A77 0IHI; iv) the left hand side is of second order d 2 X, as 
in models of thermalization and symmetry breaking in the early universe [0] • All these changes 
can be incorporated in our formalism by adaptation of the so-called Feynman rules without 
changing the combinatorial core of this paper. Hence, it is mostly for notational convenience 
that here we only consider the SPDE 

A perturbative expansion in the sense of formal power series without any control on the 
convergence of the series, as done in this work, can certainly be criticized from a mathematical 
point of view. In fact, even the existence of a solution to for the general Levy case to our best 
knowledge has not been studied though the restrictions on the nonlinearities in the continuum 
|1] that exclude polynomial interactions are related to the above mentioned UV-singularities 
and most likely can be dropped in the lattice regularized case. In this work we happily follow 
the tradition of perturbative quantum field theory (QFT) that, in spite of this, such expansions 
deliver (after suitable resummation |211 Chapters 42.5-42.7]) highly precise numerics. 

The paper consists of seven sections. In the next section we give a solution of Q in the 
sense of formal power series in the parameter A 

00 

X(t,x) = ^(-A)^(i,x). (2) 

i=o 

Here Xq is the solution of the linear equation (A = 0) and the Xj, under suitable conditions on 
the stochastic driving force 77, are determined recursively. 

In section 3 we represent this solution as a sum over rooted trees with two types of leaves, cf. 
Appendix A for the graph theoretic notions used in this article. One type of leave is standing 
for the noise rj and the other for the initial condition, cf. the following first order example for 
p = 3 

X(t, x) = x _® + x ^ + A [x-^O + X^J® + 3 + 3 X-^o] + o(A 2 ) (3) 
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In section 4 we then develop the graphical calculus for the perturbative evaluation of the cor- 
relation functions, generalizing what in physics is called the Parisi-Wu method from Gaussian 
to Levy noise. While only the well-known Feynman graphs occur in the Gaussian case, we have 
to pass on to a more general class of graphs in order to deal also with Levy type noise. We also 
prove a Linked Cluster Theorem that shows that the truncated correlation functions (cumulants 
of the process X) are given by the sum over connected Parisi-Wu graphs, only. So far no detailed 
assumptions on the stochastic driving force were needed. 

In section 5 we show how the rules for the analytic evaluation of graphs simplifies if the noise 
i] is to be taken as a white noise of Levy type. We also show that a Levy noise r) fulfills the 
requirement for the existence of the perturbative solution as specified in section 2. 

In section 6 we study the limit where t goes to infinity and give the perturbative expansion 
of equilibrium correlation functions. 

In section 7 we briefly study the question how to extract statistical information on the noise 
from an empirical correlation function, that e.g. could be obtained by X-ray analysis or sampling 
methods. To illustrate this point we solve a least square minimization problem to first order 
perturbation theory and give some indications, how to perform higher order calculations. 



In this section we give a formal solution of stochastic partial differential equation ([T]). 
We first fix some notation. Let T = R x L s = {x = (t, x) ; t G R , x G L 5 } and = [0, ?f) d 
where the opposite edges are being identified. Let B = Lg or B = V. By S(B) we denote the 
Schwartz Space of all rapidly decreasing functions on B endowed with the Schwartz topology, 
its topological dual is the Space of tempered distribution noted by S'(B). We denote by (. , .) 
the dual pairing between S(B) and S'(B). 



Now we introduce the following notation : For / G ^(r), f f(x)dx = YlxeL s ^ J*r *)dt, 
x = (t, x). For / G S{L 5 ), f /(x)dx = Ez 6 L 5 <^/( x )- Furthermore, for A C T (A C L 5 ), 
f A f(x)dx = J lA(x)f(x)dx (J A /(x)dx = f l^(x)/(x)dx) where 1a is the characteristic func- 
tion of the set A and / G S(F) (f G S(L S )). 

Let T : S(T) — > R x be the Fourier transform on S(T) 
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Formal solution of the SPDE 



i.e. 




(4) 



The inverse Fourier transform of, T ', is given by 






i.e : 




(6) 



The inverse Fourier transform of, T is given by : 




(7) 
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By definition the action, of the Lattice Laplacian A on a test function / £ S(Lg) is as follows : 

A/(x) = <r 2 [-2d/(x)+ £ /(y)]. (8) 

|x-y|=<5 

Then 

d 

^(A/)(p) = <T 2 [ - 2d + 2 £ cos(5 Pj )] Hf)(p). (9) 

i=i 

Where ej = (0, 1, 0, ...) is the canonical basis of R d and pj = p-ej. 
We define 

d 

4 m (p) = 2 ( 5- 2 (d-^cos(«5 Pj )) +m 2 . (10) 
i=i 

We introduce two convolutions product " * " and " * " respectively by 

f*9(x) = J f(x-y)g(y)dy, f,ge S(T). (11) 

and 

/ * ff (x) = J /(x - y)<?(y) dy , f,g e S(L 5 ). (12) 
Let G(i, x) be the Green function which satisfies : 



= AG(t, x) - m 2 G(t, x) + <5(rc) , 
G(t, x) = 0, f < 



(13) 



Here 5(x) is the Dirac distribution defined by 5(x) = 5(t)6(x) = S(t)S Sq x where 5(t) is the 
Dirac distribution on R and 5 Xj y = nf=i ^S- 1 x i ,S- 1 yi with <5jj the Kronecker symbol. 
Applying to the Fourier transform ,JF, it is obvious that 

HG)(E, p) = -— L . (14) 

-^ + < m (p) 

Let G*(x) be the Green function which satisfies : 

f *|(x) =A G t(x) _ m 2 ( 5 t ( x ) , 
\ G (x) = 5(x). 

Then by application of the Fourier transform T we get 

f{G t ){p) = e -^,m(P) . (16) 

By the use of the residue Theorem, one can see that G and Gj are related by the following 
equation 

G(t,x)=0(t)G t (x), (17) 

where 0(t) = 1 if t > and 0(t) = else. Let A n d be the Laplacian with cyclic boundary 
condition on [0, ?f] = Tlj. We get 
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Lemma 2.1. V N G N, 3 K N such that : 

tm 2 

e 2 

| G(t, x) |< K N : - (t, x) G ]0, oo[ x Lj. 

Proof. It suffices to prove 

sup I G{t, x) (1+ | x j 2 )^ | < K N e~ tJ ^ . (18) 



G(t, x) (1+ | x \Y = -L-(l+\ x \Ye(t) [ e^e-^l^dp 



(2 



-*^,m(P) dp 



Then 



sup \G(t, x) (1+ | x | 2 )^ | < e 

xELs 



(27T) 



/ 27T \ ' 



"2 



(1 - A u1 ) N e~ tfl l™ {p) 



dp. 



sup 



Now take 



( -T- j SUp 



" 7 te(o,oo), P enf 



1- A nf j- e 



N -tu rnip) 



V2 



tl*g _m_(p) 



(19) 



(20) 



(21) 



Definition 2.2. We define a space of all measurable and polynomially bounded functions by : 

I/I fc (*) 



Mi 



f: r 



measurable; V £; G N, 3 TV GN: 



(1+ | x 



\2\N 



dx < oo 



(22) 



Lemma 2.3. Mb is an Algebra under multiplication. 



Proof. For f,g G Mb, fg and / + ag, a G M, is a measurable function. The vector space 
structure of Mb easily follows from the triangular inequality. 
Furthermore, for k G N fixed, 3 N G N such that 



f 2k (x) 



(1+ | X | 2 ) w 
Due to the Cauchy Schwartz inequality 

\fg\ k {x 



dx < oo and 



9 {x) 



(1+ | x 



\2\N 



dx < oo 



(23) 



(1+ | x 



I2W 



dx < 



P k (x) 
(1+ I x I 2 ) 



2W 



<7 2fc (x) 
(1+ | x | 2 )^ 



dx < oo. 



(24) 
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Lemma 2.4. Let f G Mb then G * f G Mb- 
Proof. Let / G Mb- We have due to Lemma EOl 



C(N) 



yi p)jv dyi ■ ■ ■ j ( 1+ Jf_^|2)jv ^yfe 



(1+ | x | 2 ) 2fe7V 

ni/Kw) 



dx dyi ■ ■ ■ dy k 



< 4 fc G(iV) (/ ^ ' dX 



(25) 



(1+ | y | 2 )^ J J (1+ | x \ 2 ) Nk ' 

Here the right hand side is finite for iV sufficiently large. The last inequality is due to the fact 
that 

4(1+ \ x-y | 2 )(1+ | x | 2 ) 2 > (1+ | x | 2 )(1+ | y | 2 ) Vx, y G F. (26) 



Proposition 2.5. Suppose that G * rj G .Mf, a.s. 1 . Let AT(t, x) = X]^=o( — ^V-X-j(t, x) &e an 
expansion of X in the sense of formal power series in the parameter A. The perturbative solution 
of the stochastic differential equation 0) is given by 



X (t, x) = G*r ? (i ) x) + GW(x) 
X j (t,x) = G* ]T 

Bj(no,ni,...) 



P- 



no\n\\.. 



I[X?(t,x), j>l- 



(27) 



i>0 



with Xj G Mb a-s.. Here Bp(nQ, ni, ...)= {uq, n±, ...> \ nQ+n\ + ... = p; ^^ im =j — !}■ 

i>0 

Proof. We have the integrated form of the SPDE 

X(t, x) = -AG * (X p ){t, x) + G * x) + G t * /(x) (28) 
It is obvious that the solution of this equation for the zero-th order in A is given by : 

Xo(t,x) = G*»7(*,x) + G t */(x), (29) 

where G and G are the Green functions defined earlier. Now we determine the solution of higher 
order in A. We have 

oo \ 

J2(-Xyx, (t, x) + G * n(t, x) + G t */(x) 
i=o / 

/ \ 



-AG* 



ran!ni !... 

\ n ,m,...>0 u 1 i>0 

\ 710+71! + ...— p 



J\ ini X^{t, x) 



+ G*77(t, x) + G t */(x) 



1 For Levy noise with all moments this will be verified explicitly, cf. Proposition [J^D 
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By comparison of coefficients we get 

X j (t,x) = G* 77-^7- II^x) ; 3> I- (30) 

no,t»i,...>0 °" 1 "" i>0 

riQ+n 1 +...= p 
OXnQ + lXni+...=j— 1 

Xq G Mb a - s - by our assumptions. Using Lemma 12.31 and Lemma 12.41 it follows by induction 
that Xj G Mb a.s. V j G N. ■ 

3 Tree expansion of the perturbative solution 

In this section we first recall some notions of graph theory. We then recall a fundamental 
result which represents the solution - in the sense of a formal power series - of the stochastic 
differential equation (JTJ by a sum over all rooted trees. Some of the graph theoretic standard 
notions, that are used here, are collected in Appendix A and also [151117] . 

A tree T is a connected graph without cycle. We also consider graphs with different types of 
vertices. A rooted tree with root x G T, denoted by x, and two types of leaves, the leaves of 
type one are noted by ® , and those of type two by O . 

We define a partial order " < " on the set of the vertices of T, V(T) by v < w if every walk 
connecting w and x passes through v. 

We note A(i) the set of all rooted trees T with root x and two types of leaves, which has i inner 
vertices with p + 1 legs. See e.g. (J3J for the trees from A(0) and A(oo) for p = 3. 

Definition 3.1. For T G A(i), the random variable B(T, 77, x), is defined as follows: 

1) -Assign x G T to the root of the tree T. 

- Assign values y±, yi G Y to the inner vertices. 

-Assign values zi,...,Zk G T , to the leaves of type one and assign a values z^,...,z{ G L§ to 
the leaves of type two. 

2) - For every edge with two end points, e = {v, id}, assign a value G(e) = G(v — w), (v < w) 
to this edge. G is the Green function defined in (|14j) . 

3) - For the j'-th leaf multiply with rj(tj, z,-) if this leaf is of type one and multiply with f(zj) if 
this leaf is of type two. 

4) - For the I — th inner vertex multiply with the coefficient n > \(^ no - n ' y ni \ > where rii, i G N, is 
the number of rooted subtrees connected to the I — th inner vertex with i inner vertices and n' Q 
is the number of the rooted subtrees connected to this vertex with zero inner vertices and a leaf 
of type one. 

5) - Integrate with respect to the Lebesgue measure dy\ ■ ■ ■ dy p dz\ • ■ ■ dzk dz^ ■ ■ ■ dz[. 

Let T G A(j), j G N, and C(j) = {T 0) i, T 0>no T 1A , T^\, ...,Tj n . G A(i) , Y.i>o in i = 
j — 1, no + n\ + ... = p}. We construct a one to one correspondence between A(j) and C(j) 
that is given in the following way: 

We first consider the case when T G A(j) is given. Let y\ be the first inner vertex of T 
and we denote e± = {x, y±} the first edge of T. Cut the tree T at y±. Let n^, i > be the 
number of the rooted subtrees with root y\ and i inner vertices denoted by T^i, ...,Tj jrii . Then 
"-0 + ^1 + ••■ + n i + ••■ = P and J2i>o = j — 1. Conversely, let x G T and T^\, G A(i), 
be a collection of rooted trees with root y\ be given such that no + n\ + ... = p, X^i>o * n * = J — 1- 
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Attach this collection of rooted trees to the edge e = (x, y\). By definition we obtain a rooted 
tree (T, x) G A(j). For the notions of cutting and attaching we again refer to Appendix A. 
The second operation is obviously the inverse of the first one. We have established the following 
result: 

Lemma 3.2. Let T G A(j), (j G N) and 1$ i, ...,Ti n . G A(i), i G N, the application 

C : A(j) — C(i) (31) 

is a one to one correspondence and 

A : C(j) — „4(j) (32) 

is iis inverse. 

Definition 3.3. Let T € A(i) and G ^(T), we define the multiplicity of the vertex v as 
M(v) = -7T7- — I , where nii > 0, is the number of the rooted subtrees connected to v 

Q \ '''O' 1 

with i inner vertices. 

The multiplicity of the rooted tree T is M(T) = M(v). 

veV(T) 

For / G S(L 5 ) we put F f (x) = <5 (i)/(x) , x = (t, x) G L 5 

Lemma 3.4. Zei T G -4(j) and fo G V(T) the vertex connected to the root x, then the following 
result hold : 

B(T,x, v ) = M(v )G*[ H B(S,., V )](x). (33) 

sgC(t) 

Where C(T) is the cut (cf. Appendix A) of the tree T. 
Proof. 

(g) dv 

e£E(T) ZieLi(T) leL 2 (T) v£V(T)\{x} 



B(T, x, n) = M(T) I H G(e) H r,(h) ]J F f (l) 

eS-E(T) ZiGLi(T) 1&L 2 (T) 

= mm j g(x - v ) n [ m (s) j n G ( e ) n ^ 



seC(T) " 17 ees(s) heLi(S) 

x II F /W ® dv]dv = M(v )G*[ ., 77)](a?). (34) 

ZeL 2 (S) t'e^(5)\{i>o} SGC(T) 

■ 

Using Proposition 12.51 and Definition 13. II we get the following theorem: 

Theorem 3.5. The solution of the stochastic differential equation (1) in the sense of formal 
power series is given by a sum over all rooted trees T G A(j) that are evaluated according to the 
ruled fixed in Definition VJ.l\ i.e. 

Xj(x, V)= Yl B ( T > 00 > 7 ?) ' 00 e r ( 35 ) 
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Proof. We prove the assertion by induction. For j = 0, we have Xq(x, rf) = G*rj(x) + G*Ff(x), 
which just is the sum of the evaluation of the two trees in *4(0) . 
Suppose that (|3*5|) is true for all i < j. By proposition (2.2) we have 



X j (x,r,) = G* Y, 

Bp (no, rai,...) 



no!ni!.. 



i>0 



(36) 



Then using the induction hypothesis and Lemmas 13.21 and 13.41 we get 



Xj(x, V ) = G * 

B J p (n ,ni,...) 



Pi 



n 



nn\n\\... -*■ 



E B ( T --^) 



(X) 



E E 

Bj(n„, ni) ...) T l>1 ,...,T l ^eA(i) 

i=0,l,2,... 



pi 



rto!n.i!. 



G*nn%.-.i) ^ 



i>0«=l 



E E s (^( I b,i,---,7b )ri0 ,Ti il ,...),a;,?7) 



Sj(n ,ni,...) T i,i--' T i,n i e-4(0 

i=0,l,2,... 

E B(T,x,f7). 



(37) 



4 A graphical calculus for the correlation functions 

In the following we define a class of graphs that we call Parisi-Wu graphs. We then give a rule 
which assigns a value to such a graph. We prove a theorem which represents the moments of the 
solution given in the previous section in terms of a sum over all generalized Parisi-Wu graphs. 
We finally prove a linked cluster theorem for Parisi-Wu graphs. 

Definition 4.1. A generalized Parisi-Wu graph of m-th order with n exterior vertices x±, ...,x n 
of fertility one, is a graph that contains n rooted trees T\ € A(ji), ...,T n £ A(j n ) as subgraphs, 
such that any leave of type one is connected to exactly one vertex of a new type, called inner 
empty, of arbitrary fertility. We denote such a graph by Q and the set of all generalized Parisi- 
Wu graphs of m-th order and n roots by P(m, n). 
If we denote by Q = Q(Q) the set of the empty vertices then 

n n 

V(G) = (j V{T k ) U Q, E{Q) = (J E(T k ) U E(T k , Q) (38) 

k=l k=l 

where E(T k , Q) is a set of edges e = {v, q}, v € Li(T k ), q 6 Q such that p(v) = 2, Vi> € Li(T k ). 
By definition, inner vertices are distinguishable and have non distinguishable legs whereas empty 
vertices are non distinguishable and have non distinguishable legs 2 . 



2 The notion of distinguishable and non distinguishable legs are explained in Appendix A. 
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The name generalized Parisi-Wu graph has been chosen to distinguish the graphs used here 
from Gaussian Parisi-Wu graphs implicitly defined in |18j . which, in a topological sense, coincide 
with the ordinary Feynman graphs. 



Proposition 4.2. For j\, j n G N, let T\ G A(ji), T n G A(j n ) n rooted trees with roots 
Xi, ...,x n G r, then the following result holds in the sense of formal power series 

In \ oo 

(T[X{x h rj)\ = (Bi(T 1 ,x 1 ,ri)'--B n (T n ,x n ,ri)) (39) 

\i=l / m=0 j'i,-,Jn>0 

3\ + ---3n—m 



Proof. 



mxix^r,)) = (j2\iX j (x 1 ,r,)--- 1 £\iX j (x n ,r,)) 

\i=l I \j=0 j=0 I 

oo 

= E A ™ E (Xnixuri.-.x^T,)) 

m=0 ii.— jn>0 

ix+..-+jn=m 

oo 

= E ( B l( T l> Xl,V)---Bn(T n , Xn , T})) (40) 

m=0 ii,.-.jn>0 
31 + ■ ■ ■+jn=m 
T^AUl), -,T„£AUn) 

The last equality is an immediate consequence of Theorem 13.51 

■ 

Let J C N be a finite set. The collection of all partition of J is denoted by T>(J). A partition is 
a decomposition of J into disjoint, nonempty subsets, i.e 

/ G V(J) 3 A; G N, / = {/ l5 ...,!*}, I, n /, = V 1 < j < I < k , U? =1 J, = J. (41) 

Definition 4.3. Let z\, 4 £ T, 1 a partition of the set {1, n}, I G £>({1, n}), 

I = {Ji, Ifc} the truncated moments functions (77(2:1) • • • r](z n )) T are recursively defined by : 

n k 

<n »?(*)>= e n^) T ( 42 ) 

1=1 76X>({l,...,n}) /=1 

where <Ij> T = <I[ ^)> T - 

Definition 4.4. For Q G P(m, n) and xi, . . . ,x n G T we define the number V\^\{x\^ ■■■,x n ) as 
follows: 

1- Assign values x±, ...,x n G T to the roots of the trees Ti, T n . 
-Assign values y\, y m G T, to the inner vertices. 

-Assign values z±, G T , to the leaves of type one and assign values z' 1; z' t G L$ to the 
leaves of type two. 
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2- For every edge in a tree T r , r = 1, . . . , n with two end points ,e = {v, w} , v < w, assign value 
G(e) to this edge. 

3- For each inner empty vertices with I legs connected to leaves with arguments q\ , . . . , q\ multiply 
with 

(ri(qi)---ri(qi)) T . (43) 

4- For the j-th leaf of type two multiply with /(z^). 

5- For each vertex v £ V{G) that is an inner vertex of a tree multiply with the corresponding 
multiplicity M(v), (see DefESJ). 

6- Integrate with respect to the Lebesgue measure dy± ■ ■ ■ dy m dz\ ■ ■ ■ dz^ dz[ ■ ■ ■ dz\. 

Let m, n £ N, I £ 2>(U"=i Li(Tj)), I = {h, ...,4} and Ti £ A(ji), T n £ A(j n ), such 
that ji + ...+jn = Tn. We construct a one to one correspondence between pairs (Ti, T 2 , T n ; I) 
and the generalized Parisi-Wu graph Q £ P(m, n) in the following way: 

Let the Parisi-Wu graph Q is given and we construct (Ti, T n ; I). By definition (/ contains 
n rooted trees, T\ £ wA(ji), T„ £ A(j n ) such that ji + ... + j n = m. Suppose that there are k 
empty vertices in the graph Q. Give an arbitrary number I = 1, k to each empty vertex. For 
the Z-th empty vertex let // be the subset in (Jj=i -^i(^j) corresponding to the leaves connected 
to this vertex. Then / = ...,Ik} £ f (U^=i ^i(^i))- 

Conversely, let Ti, T n , I be given such that Ti £ »4(ji), T n £ *4.(j n ) and / £ T>((J™ =1 ^i(^j))> 
I = {ii, ifc}. Let Q = {!;•••>&} be the set of the empty vertices of £/. Let V((?) = 
5^" =i y(Tj) U Q be the set of all vertices of Q and E(Tj, Q) be the set of all edges which 
connect the j — th rooted tree Tj, (1 < j < n) with the empty vertices in Q. Then E(Tj, Q) = 
{{v, q) : »e Li(Tj), q £ Q ; u £ /J. Let = E"=i ^( T j) u Q) be the set of a11 ed g e 

of Q. 

We have established the following result : 

Lemma 4.5. Let m, n £ N and M(j) = {Ti £ AO'i), r n £ «4(i„) ; I £ X>(£" =1 T(Tj)), I = 
Ifc}}. T7ien i/ie mapping N constructed above 

N : P(m, n) — ► M(j) (44) 

is a one to one correspondence. 

Theorem 4.6. Te£ {X t , t > 0}, 6e the solution of the stochastic differential equation (1), and 
Ti,...,T n , n rooted trees with roots x\,...,x n . Then the moments of X(x) are given by a sum 
over all generalized Parisi- Wu graphs Q € P(m, n) of m — th order that are evaluated according 
to the ruled fixed in definition *- e '■ 

In \ 00 

\i=l I m=0 5eP(m,n) 

Proof. Let 



n n n n 

Li(Q) = [jL 1 (T ji ), L 2 (G) = \jL 2 (T ji ),V(G) = (JV(T ji )\({x 1 ,...,x n }), E(Q) = \jE(T ji ) 

i=i i=i i=i i=i 
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The multiplicity of the graph Q is M(Q) = M(i 

v e V{g) 



we have 



Y[X( Xi ,ri)\ = ^A™ Yl B^Tuxu r)) ■ ■ ■ B n (T n , x n , rj) 

U=l / m=0 \ n,...,jn>Q / 

31+—+jn=m 
T^Aiji), ...,T n eA(j„) 

oo 

= e Am e / n 

m=0 7i,-,Jn>o ^ ee£(S) 

Tie^(ii), .-,T n e-4(jn) 

x n F fQ)( n ) ®«6V(o)\g(g) ^ 

ie£ 2 (a) \Jieii(s) / 

oo „ 

= e Am e m (£) / n ^ 

m=0 ii,-..,jn>0 J e eE(g) 

il+.-.+in—m 

Tie^Ui), ...,r„e-4(j'n) 



x II F f0) E IIvII^ 1 )) ®vev(G)\QM dv 
leL 2 (g) ieT>{Lx(g)) 1=1 \heli I 

/={/i,-..,/ fc } 

oo . 

= e Am e / n c ( e ) 

m=0 n,-,Jn>0 ^ e&E(Q) 

il + ---jn—m 

I6Z>(i!(S)) 
J={J 1 ,...,/ fe } 

x n F /wn(n*)) ( 46 ) 

«gL(6) 1=1 \heii I 
Now we apply Lemma 14.51 and Definition 14.41 to conclude. 

■ 

We denote the collections of the connected generalized Parisi-Wu graph of m-th order with 
n roots by P c (m, n). 

We construct a one to one correspondence between P(m, n) and Af(m, n) = {I G V({l...n}), I = 
{h, 4} ; G\ G P c (m, pi), G P c (m k , p k ), m = mi + ... + m k , n = pi + ... +Pfc} in the 

following way: 

We first consider the case when a generalized Parisi-Wu graph Q G P(m, n) is given. By 
definition £/ contains n rooted trees T\ G *4(ji), T n G x4(jn), with roots ;ci,...,x n and an 
arbitrary number of empty vertices. For 1 < q < j < n we say that the graph connects q 
and j, in notation q ~r j, if x g and Xj are connected in the generalized Parisi-Wu graph. 
Obviously ~# is an equivalence relation on {l,...,n}. Let I = {Ii, ...,/&} be the equivalence 
classes of ~# then I G T>({l...n}). For 1 < I < k let Q\ be the connected generalized Parisi- 
Wu graph with m\ inner vertices and pi = j{Jj roots, then m = m\ + ... + m k and n = p\ + 
... Conversely let 7 G £>({l...n}), I = {h, 4}, Gi G P c (m, pi), G P c (m k , p k ), 

m = mi + ... + m k . 
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Obviously the connected generalized Parisi-Wu graphs Gi,---,Gk determine a generalized 
Parisi-Wu graph Q with m = mi + ... + inner vertices and n = p\ + ... + Pk roots by taking 
the disjoint union. 

This second operation is the inverse of the first one. We have established the following result: 
Lemma 4.7. The map 

M : P(m, n) — ► J\f(m, n) (47) 

is a one to one correspondence. 

The following is known as linked cluster theorem and gives the truncated moments of X(x) 
as a sum over the connected Parisi-Wu graphs, only. 

Theorem 4.8. Let P c (m, n) be the set of connected generalized Parisi-Wu graphs, then the 
truncated moments [X(x\, rj) ■ ■ ■ X(x n , rj) T are given by 

oo 

(X( Xl ,r,)-.-X(x n , V ) T = ^A m £ V[g](x u x n ). (48) 

m=0 geP c (m,n) 

Proof. By Theorem (|4.(ij) we have 

00 

(X(x 1 , V )--.X(x n , V ) = J2 Xm E nO}(x u ...,x n ). (49) 

m=0 £eP(m, n) 

In the other hand by definition of the truncated moments 

k 

(x( X1 , V )-..x(x n ,r ] ))= yi n<n x ^' ? ?)> T ( 5 °) 

I6T>({l,...,n}) 1=1 jel t 

I={I U . ..,-r fc } 

To prove the theorem it suffices to prove that in equation (|5U[) we can replace the truncated 
moments by the right hand side of ([48 [I. i.e. 

k I oo \ 

(x( Xl , V ) ■ ■ ■ x( Xn , V ))= II E xmi E ^[fti^i ; 3 G ft ■ ( 51 ) 

;ec({i,..,n)) 1=1 \m 1= ftePcC^JA) / 

/ ={ / l I k! 

Thus we have to prove that the right hand side of 1)51(1 is equal to the right hand side (r.h.s) of 
l|39"|). We have for the r.h.s of (|5*T1) 

oo oo 

E E Ami "- E P[*;iG4) 

m 1= m fc=° 

i={i lt ...,i k } SiePWctmi.m) S fc eP c (m fc , Pfc ) 
oo 

E E A -•""7n;,T /,!••• V[Q k ]{xj ; j € h) 

I ET> ({1,..., n}) mi,...,m^=0 
/={/!, ...,/ fe } eiePc(roi.Pi) 

e fe 6P c (m fc ,p fc ) 
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oo 



£ X m 1+ ... +mk 



rn\,...,m k =0 Jez>({l, ...,»}) 

l={l lt ...,l k } 

Si ePc (mi, Pi) 



OO 

E 

m=0 Jex>({l, ...,«}) 



A m £ ^i]fe;i6A)--»j;ie4) (52) 



/={/!,.. .,Ife} 
Gl6Pc(mi,pi) 

G k ePc(™ k .p k ) 

m± ,...,m k — 0,1, .. . 



Now application of Lemma (|4.7|) concludes the argument. 



5 Levy noise 

In this section we define first a white noise measure with Levy characteristic tp, then we recall a 
theorem which give the truncated moments of the noise which permit us to derive a simplification 
of the Feynman rules (Definition I4.4JI . 

Let v be an infinitely divisible probability distribution. By Levy-Khinchine theorem (see e.g. 
0) we know that the Fourier transform (or characteristic function) of v, denoted by C u , satisfies 
the following formula 

C„(t) = f e ist du{s)=e^ t \ teR 
Jr. 

where ip : R — ► C is a continuous function, called the Levy characteristic of is, which is uniquely 
represented as follows 

ijj(t)=iat- — + [ ( e ist - 1 - dM(s), Vfe M. (53) 

2 Jr\{o}\ 1 + s 2 J 

where a £ R, a 2 > and M is a Levy measure on R \ {0}, i.e. f R ^ Q -^ mm(l, s 2 )dM(s) < oo. 
Under the additional assumption that all moments of M exist, i.e. Jm/Q} s 2n dM(s) < oo, 
n € No, one can reparametrize this representation setting z = J^^ ydM(s), r = M/z, a = 
^ ~~ Jk\{o} J^?dM(s) and one obtains 

gH 2 r 

if){t)=iat + z (e lst - l) dr(s), V t G R. (54) 

2 Jr\{o} 

The first term in Q54|) is called deterministic, the second one Gaussian term and the third one 
Poisson term. 

We note by B the cr-algebra generated by the cylinder sets of S'(T). Then S'(T) is a measurable 
space. 

We define a characteristic functional on S(T), as a functional C : S(T) — ► C such that : 



1. C is continuous on S(T); 



How to determine the law of the noise driving a SPDE 



15 



2. C is positive-definite; 

3. C(0) = 1. 

By the well-known Bochner-Minlos theorem (see 1121 ) there exists a one to one correspondence 
between characteristic functional C and probability measures [i on (S'(T), B) given by the 
following relation 

C(f)= [ e^d^), f G S(T). (55) 

JS'(T) 

We set r± = {x G T : x = (t,x), ±t > 0}. 

Theorem 5.1. Let ip be a Levy characteristic given by the representation \5ty then there exist 
an unique probability measure P^ on (S'(T), B) such that 

C v (f) = [ e ir '^dP^{i) = exp ( [ i>(f(x))dx) , x = (t,x). (56) 
Js'(r) \Jr + J 

Where rj(f, = (/, C), f e S(T). 

Proof. The right-hand side of (|5(ij) is a characteristic functional on S(T) (see eg. |12j). The 
result thus holds by using the Bochner-Minlos Theorem. 

■ 

Definition 5.2. We call P^ in Theorem 15. II a generalized white noise measure with Levy char- 
acteristic ip and (S"(r), B, P^) the generalized white noise probability space associated with i/j. 
The associated coordinate process 

rj : S(T) x (S'(T), B, P*) — . C, r,(f)(u) = u(f) V/ G 5(T), w G 5'(r) 

is called a Levy noise. 

The following result on Levy noise is essential for our perturbative approach 

Proposition 5.3. There exists a version of G * n such that G * rj G Mb Pip &- s - ■ 

Proof. The proof requires several steps. 

i) By (|54|). the Levy noise can be decomposed into independent parts n = if + rf + rj p , 
where rf, rf and rf are respectively the deterministic, Gaussian and Poisson noise given by 
their characteristic function in theorem 1)5.11) . It is sufficient to deal with these parts separately. 
The statement holds trivially for the deterministic part. 

Let us next focus on the Gaussian part. We get the following two results that together show 
G*t] 9 G M b a.s. . 

ii) G * rj 9 (t, x), for x G L$ fixed, has a continuous extension in t. In particular, G * rf is a 
measurable function on T (a.s.). 

By the Linked Cluster theorem we have 

(| G*rj 9 (t,x) -G*rj 9 (s,x) | 4 ) = 3((| G * rj 9 (t, x) — G * rj 9 (s, x) | 2 )) 2 . (57) 
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If we assume that t > s, then it is easy to see, using (r] 9 (t, x)r/ 9 (s, y)) = a 2 9(t)9(s)5(t — s)<5(x— y), 
that 

(| G*rj 9 (t,x)-G*ri 9 (s,x) | 2 } 

= ((G * v 9 (t, x)) 2 + (G * rj 9 (s, x)) 2 - 2(G * ^(t, x))(G * ^(s, x))} 

2 



+ 



a 



9{t- t')G t - t > (x - x')r] 9 {t', x')dx'dt' 

d(s - i')G s _ t /(x - x')rf(t',x')dx.'dt r 

9(t - t')9{s - J)G t -v (x - x')G s _ s ,(x _ x ") 

x ?f (i', xV(s', x")dxdx"dt'ds' 



9{t - t')Gt t ,(x ~ x') + 0(* - t')G 2 s _ t ,{x ~ x) 



/o JO 



29(s - t')G t -t'(x ~ x')G s _ t /(x - x) Weft' 



a 2 / / e(s-t') G t _ t ,(x') -G s _ t .(x') dx'dt' + / G 2 _ t ,{x')dx'dt' (58) 




Moreover by the expression of Gt given in (|16j) we get : 

2 C 2 (t-S) 2 



(l + \x\ 2 ) d [G t (x)-G s (x 



< 



(27T) 



2</ 



and 



(l + |x|YG 2 _ t ,(x 



-2d 



< C 



(27T) 



2J ' 



(59) 



Here c = sup 4>0 J" n d(l — A) d e ^ m ^dp < oo. Inserting these two estimates in (|58j) . we get for 



X > sufficiently large 



G*rj 9 (t,x) -G*rj 9 (s,x) | 4 ) < K \ t - s 



(60) 



Application of Kolmogorov's extension theorem now proves the assertion. 

ii) G * r] 9 (t, x) is polynomially bounded (a.s.). 
We prove the stronger statement that the expectation of J dx is finite. We have for 

k 6 N even 



(G*7?g(x)) fc 
(1 + |x| 2 )^ 



.4 



G*rj 9 (x) | 2 ) 



(1+ | x 



\2\N 



dx 



< Aa z / G z (x)dx 



dx 



(1+ | x l 2 )^ 



< oo. 



(61) 



The last inequality hold for an arbitrary N > d/2. Here A is the number of pairings of k objects. 
It is left over to deal with the poisson part. First we recall the path properties of rf following 

PI- 
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iii) rf is a signed measure with locally discrete support (a.s.). 
Let A n C r + be a monotone sequence of compact sets s.t A n j T + as n — > oo and Aq = 0. 
For n € N let D n = A n \ A n _i and we denote the (Lebesgue) volume of D n by | D n |. Let 
f]n = J2j=i SJ^Xf be a random field, where 5 X is the Dirac measure of mass one in x and N* is 
a Poisson random variable with intensity z \ D n | i.e 

P( N * = l) = e -*IA,| (*l Dn\) 1 . l € Nq _ 

For n € N, {X™}j g j^ is a family of i.i.d. T-valued random variables distributed uniformly on D n . 
{■5?}i,neN is a family of real valued random variables with law given by r. All these random 
variables are independent of each other. The characteristic functional of f% is given by 

,ir&(f)\ = / e i^ksp{Xf)\ 

-z\ Dn \ j2 i^t (f f j>m dr (s" dx X 



i! \JD n JR\{0} \ D n\J 

exp \z [ I (e is ^ - l)dr(s)dx 1 = C vP (f), V/ G 5(r),su PP / C D n . 



(62) 



Hence rf ', when restricted to D n , coincides in law with rfi. By the uniqueness statement of 
the Bochner-Minlos Theorem and the fact that by construction f% takes values in the locally 
discrete signed measures (a.s.), it follows (cf. [2] for the details) that also rf with probability 
one is a locally discrete signed measure. 

iv) For / € L 1 (r,dx) bounded and non-negative, rf is /-finite (a.s.), i.e. J r f d\r] p \ < oo. 
Here \rf\ = rf. + rf. is the modulus of the signed measure rj, cf. [T3"] . 
This assertion can be seen from the following calculation 



L 



fd\rf\) = £;// fd\ff n 



n=0 
oo I N% 

E(E I s ?\fW 

n=0 \ j=l 



^ ^ C* - !) ! I I Jr\{0} 

y^ e -*|D n | e *|D n | / f( x ) dx \ s \dr{s) 

„_n ■/£>„, ./RUO} 



2 

n=0 



'M\{0} 

s | dr(s) I f(x)dx <oo. (63) 



M\{0} 

v) G * rf 6 -^ 1 (r, ft da?) a.s., where g e = (1 + |a?| 2 ) (c?/2+e) _ j n particular, G * rf is 
measurable. 

Let D l n = A n \ Ai for n > I, we denote the restriction of the noise rf to an open set A C T 
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by 7/j\ Clearly, G * ryj^ € -^(T, g e dx) (a.s.) since G is in L x (r, dx) and supp^ is finite 
(a.s.). The following estimate shows that G * 77^ forms a Cauchy sequence in L 1 (r, g t dx). 
With || . || €j i the L 1 — norm on that space, we get 

sup||G*??f -G*rff \\ e 1 = sup || G*r/f || e 1 

n>l lAn |A < n>Z H 

< sup / I G I *g e d I 7/f I 

n>lJT ,D n 

= / I G I * g e d I r/ p | — > as I — ► 00 (64) 

since ?7 P is |G | * ^-finite (see iv) and note that L 1 (r, dx) is closed under convolutions). Also, 
lim (G*rf An , /)= lim G* /) = {rf , G* /) = (G*rf , /) V/ € S, (65) 

n — > 00 " n > 00 n 



and by the fact that convergence in L 1 ^, g e dx) implies convergence in S', we get that G * if 
coincide with the limit of G * rf^ in the Banach space L l (T, g e dx). 

vi) G * r] p (t, x) is polynomially bounded. 

({G*T]P(x)) n ) f ^ -X ((G*f] p (x))^'' rj 



J (1+ I - \ 2 ) N ' J ,J? _ 11 (1 



__|_ I x \2\N 

«}) 1=1 

I={/l,...,/fe} 



<7.r 



/6B({l r ..,«}) J /=1 V 11/ 

S C/— ^ F w<». (06) 



The last inequality hold for an arbitrary N > d/2. Here 

n 

C = b n max sup ' J | ((G * ?] p (t, x)) n; ; 
n i-T,i=i n i= n (*,x)er l=1 



T 



with 6 n the n-th Bell number, i.e. the number of partitions of n objects, and the supremum of 
the truncated expectation values is finite by virtue of Theorem 15 . 41 b elow and Lemma 12. II ■ 

Theorem 5.4. The truncated moment functions of the Levy noise r/ are given by the following 
formula 

«*,)■ f *(.-«)-. (67) 

where 



r-, 



y ' dt n 1 



= 5 n . 1 a + 6 n ^o- 2 + z s n dr(s) (68) 

JR\{0} 

<5 n „/ fremg £/ie Kronecker symbol. 
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x — <— o<K>>- x x — < — O — > x 



Figure 1: Construction of a generalized Feynman graph where every leave of type one 
together with the two edges connected to it is replaced by one edge 

Proof. After application of the linked cluster theorem to (|56|) . the formula follows by a straight 
forward calculation, cf. for details. ■ 



The particular form of (|67|) makes it trivial to carry out the integrals over the leaves of type 
one of the Parisi-Wu graph. This leads to the following simplification of Feynman rules: 

Theorem 5.5. Suppose that the noise r] is of Levy type. Note that in a Parisi-Wu graph, each 
leave of type one has exactly two legs. For a given Parisi-Wu graph Q € P(m,n), one can thus 
get an equivalent graph Q' where every leave of type one together with the two edges connected 
to it is replaced by one edge (cf. Figure ^)). 

Let P' c {m,n) be the collection of all connected graphs obtained in this way from P c (m,n). 
The Feynman rules V[Q'] for the graph Q' then simplify in the following way: 

1- Associate the values xj, j = 1, ...,n, to the roots, an integration variable y to each inner 
vertex and an integration variable z to each leave of type 2. 

2- For each inner vertex v of a tree Tj, multiply with a multiplicity factor M{v) and for each 
inner empty vertex with I legs multiply with a factor q . 

3- For each leave of type two multiply with a factor Ff(z). 

4- For each edge multiply with a propagator G(e). 

5- Then integrate over all inner vertices and all leaves of type two. 
In this way, one obtains in the sense of formal power series 

oo 

(X( Xl )...X(x n )) T = ^(-A) m nG']{x 1 ,...,x n ) (69) 

m=0 g'ePt(m,n) 

6 Equilibrium correlation function 

In this short section we determine the limit of the correlation functions when all time arguments 
go to infinity simultaneously. 

Definition 6.1. Let P[ c (m, n) be the collection of graphs in P' c {m, n) that do not posses a leave 
of type 2. For Q £ P[ c (n,m) the value Voo[^](xi, x n ) is obtained the following Feynman 
rules: 

1- Assign values (0,xi), (0, x n ) G L$ to the roots of the trees T\, T n . 
-Assign values y±, ...,y m € T, to the inner vertices of the tree. 

-Assign values zi, Zf. € T to the inner empty vertices. 

2- For every vertex, v in a tree, multiply with the multiplicity coefficient M(v) and for each 
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inner empty vertex witl I legs by q. 

3- For every edge with two end points e = {v, w}, (v < w), assign a value G(e) . 

4- Integrate over T_ with respect to the Lebesgue measure dy\ ■ ■ ■ dy m dz\ ■ ■ ■ dzk- 

Lemma 6.2. Let Q G P' c (m, n) and let I{Q) = Y\ e &E(g) G ( e )> then ViV G N, 3 if = K(N, m) 
such that 

i m i < 7— • ( 7 °) 

( 1 + max \v — w \ ) 

v,w€V(g) 

Proof. Let v',v" G V(Q) such that I v' — v" I 2 = max I v — w I 2 . As Q is connected, 

v,wev(g)\{v } 

there exist a walk W from v' to v" and let q < m + 1 be the number of steps from v' to v" . Let 
e±, e2, e g be the edges of the walk W. Then 

= n c ( e ) x n 

e&W eeE{Q)\W 



In W, there must be at least one of the q steps that is > ^ ■ . By the use of Lemma 12. II 
we find 



K(N,m) , , 



Theorem 6.3. Zei Q G W c (rn, n), then the perturbation series for the truncated moments 
converges in the sense of formal power series, when t goes to infinity to : 

(X 00 (x 1 )---X 00 (x n )) T = lim (X(xi) ■ ■ ■ X(x n )) T (73) 

t > oo 

where 

oo 

(Xoo(xi) • • • I„(x„)) T = Xm E VooPKxi, ...,x n ). (74) 

m=0 GeP[ c (m,n) 

Proof. Let P2 c ( m i n ) = ^°c( m > n ) \ A' c( m > n ))- We have 



(X(x!)-X(x n )) T = V[g}( Xl ,...,x n ) + Y^ m E V[0](x l5 ...,x n ) 

m=0 GeP{ c (m,n) rn=0 GePt, c (m,n) 

As we take the limit i — ► 00 in the sense of formal power series, it suffices to show 



(75) 



lim V[e]((t ) x 1 ),.,(t ) x n ))=V m [e](xi,...x„),VgeP( c (m,rj) (76) 

t — >oo 

and limt—Kx, V[G]((t,x.i), (t,x n )) = 0, V£? G P^d 171 ^ n ) 
Let first G <E P{ c (m, n), we have 

V[0]((*,X!), ...,(*, x n )) = M(0) / TT G(e) <fo (77) 

1/1 (o,t) eGE(e) «ey(e)\{(t ) xi),...,(t,x„)} 
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where I\ a M = { x = (^ x ) S T,a < t < b}. Now we transform s — > s — t for v = (s, v) E V(£7) 
and let t go to infinity. 

Let A e = {e = {n, «;}, u = (t v , v), w = (t w , w), (t> < u;)}, we get 

lim V[g]((i,xi),...,(t,x n )) = Af(0) lim / ]~[ G ( e ) ® dv 

t > oo t— >oo /™ — 

1/1 (-*,o) A e vev(g)\({ Xl ,...,x n }) Xl= ( 0jXl ) 
= V 00 [g](x 1 ,...,x n ) (78) 

where the convergence in the last step is due to Lemma 16.21 

Let now Q £ P^^m, n). Then, again by Lemma 16.21 1(G) is rapidly decreasing in the 
difference of the time argument t of the external vertices and the time argument of the leave 
of type two. Consequently, the integral over 1(G) gives no contribution in the limit t — > oo. ■ 



7 Determination of the law of the noise 

The main subject of this section is to find statistical information on the law of the noise r\ driving 
the SPDE (Q) from empirical data. One approach to achieve this is the least square method, 
i.e. to solve the minimization problem for the two point function of the stationary distribution 

[ Q(ci, c 2 ,...) = / | F th (x, ci,c 2 ,...) -F em (x) | 2 (he, 
(^=0,jGN, (79) 

where F em is the empirical 3 correlation function. Note that by Lemma (|6,2j) the function 
Fthfe) = (X oo (0)X oo ('X-)) T to any order of perturbation theory is rapidly decreasing in x and 
thus there should be no problem with the convergence of the integral in (|79j) - if the modeling 
is not completely wrong, F em should also be of fast decay. 

In the remainder of the section, we give the solution to (|79|) to the first order in perturbation 
theory for p = 3. To simplify the calculation we note first 

Lemma 7.1. Let the measure r be symmetric, i.e r(A) = r(—A) V A G B. Then one can omit 
all such generalized Parisi-Wu graphs from the perturbation series that have an empty vertex 
with an odd number of legs. 

We assume this symmetry of r in the following. On inspection of the fist order solution, we 
obtain by the use of Theorem (j-5,5|) : 

F th (x, c 2 ,c 4 ) = c 2 P 1 (x) + Ac 4j P 2 (x) + (A 2 ) (80) 

where Pj(x.) = Voo[Gj](0, x)/c2j, j = 1,2 and Gj the first/second graph in the above first order 

expansion (cf. Figure EJ. Note that Pj does not depend on c 2 ,C4 anymore. 

Apparently, Q only depends on C2 and C4. We have to solve the equations dQ(c2, c 4 ) = , £ = 2, 4 

Moreover, 

Q(c 2 , C4) = a c\ + (3 c\ + 27 c 2 c 4 - 2 c 2 a - 2 c 4 b + c (81) 

3 In material sciences, the empirical correlation function F em (t, x, c\, C2, ...) can be measured by X-ray 
spectroscopy or sampling methods, cf. e.g. [51ll4|. 
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Fi h ,(x, c 2 ,c 4 ) = x o x + x «Q0 x+ °( A2 ) 



Figure 2: Expression of the correlation function in the first order expansion. 



where a = J Pfdx, = A 2 / PfcZx, 7 = A /PiP 2 ^x, a = $ PiF em cht, b = Xf P 2 F em <ix and 
c = f F% m dx.. We now write down the equations =0,1 = 2, 4. and solve for c 2 and C4. 
One then obtains the following first order approximation of c 2 and C4: 

a 13 — 7 6 . , 



C4 = 7-^5 2 (83) 



and 

06 — 70 

a (3 — 7 2 

c 2 gives a measure for strength of the fluctuations of r\. Let the kurtosis K given by the following 
expression: 

K = °4 (84) 

If K = , then there is no jump, the stochastic dynamic is purely diffusive. If < \K\ « 1 , 
there are some jumps but the stochastic dynamic is predominantly diffusive. If \K\ » 1 , then 
the stochastic dynamic is predominantly ruled by jumps. 



8 Appendix A 



For the convenience of the reader we collect some graph-theoretic notions which have been used 
in this work. 

Let V be a finit set and 

E = {e = {v, w} ; v, w G V}. 

The elements of V are called vertices, which are of different types in our case we have the roots 
x , the inner vertices • , the empty vertices o and the leaves of type one ®. Such vertices are 
labeled by arguments in T. We have another type of vertex, called the leave of type two, O , 
which are in L$- 

The elements of E are called edges , i.e lines connecting exactly two vertices, and we say that 
an edge e = {v, w} joins v and w. Thus we can define a graph 

Q = Q{V) 

with the vertex set V as a family of pairings 



E(Q) C {e = {a, b}, a, b £ V}. 
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For a given graph Q we note E{Q) (resp. V(G) = V) the set of edges (resp. vertices ) of Q. 
A graph H is called a subgraph of the graph Q and we write H C G when the vertex set V(H) 
of H is contained in the vertex set V{G) of Q and all edge of H are edge in Q. 
For vi, v n € V(G) we define a v\ — v n walk on Q as a sequence of vertices and edges 
W = (v\, ei, U2, e2, ...,u n ) such that 

if fi = u n we say that is a closed walk on £/, this closed walk is said to be a cycle when 
e { / e 3 ; V i,j G {l,...,n- 1}. 
degree of a vertex. 

Let <7 be a graph and u € V(G), the degree of u, noted by p(v) is defined as : 

p(v) = He e E(G) : e = {v, w}, w € V} + fj{e € E(G) : e = {v}} 
connected graph. 

A graph G is called connected if there is a v — w walk for all v, w G V(G), otherwise G is 

disconnected. 

rooted tree. 

a rooted tree is a pair (T, u) such that T is a tree and v a vertex of T with p(v) = 1. v is called 

the root of the tree T. 

leaf. 

Let (T, «) be a rooted tree, then any vertex w of T such that p(w) = 1, w ^= v, is called a leaf 
of T. 

We note E(T), V(T) and L(T) respectively the sets of edges, vertices and leaves of the rooted 
tree T. 

rooted subtree. 

A rooted subtree of the rooted tree (T, x) is a pair (5, v) where S is a subgraph of T and 
v £ V(T) such that Vm S V(5) we have v < w and p(v) = 1. 

A cut, C(T) of the rooted tree T in the vertex u, directly connected to the root by an edge 
e = {x, v}, is the uniquely defined collections {Ti,...,T n } of rooted subtrees of T with root v 

n 

such that 1J L(Tj) = L(T). 
i=i 

Attachment of rooted trees. 

Let x 6 r and (Ti, y)...(T n , y) be a collections of rooted trees with V{Ti) and E{T{) design 
respectively the set of vertices and edges of the i— th rooted tree. We define an attached rooted 
tree (T, x) by : 

n n 

V(T) = V(Ti) U {x} and -E(T) = E{Ti) U {x, y} respectively the set of vertices and edges 

i=l i=l 

of T. 

Rooted trees with two types of leaves. 

Let (T, x) be a rooted tree and L(T) be the set of the leaves of T, (T, x) is said to be a rooted 

tree with two types of leaves if and only if : 

L(T) = Li(T) U L 2 (T) and I £ L{T) \L x {T)^l £ L 2 {T). 

We said that L\{T) is the set of the leaves of type one of T and L 2 {T) is the set of the leaves of 
type two of T. 
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(T) 

In this graph (T) the inner vertices are distinguishable and have non distinguishable legs : 
Let x be the root of the tree T and we consider the edges e\ = {y±, z±} and e 2 = {yi, z 2 }, the 
value of the tree T [cf. Def.(3.1)], is the same when we permute the edges i.e when e\ = {y\, z 2 } 
and e 2 = {yi, z\}. 

In this Parisi-Wu graph (£/),with one inner vertex y\ and two roots x±, X2, the empty vertex 
have non distinguishable legs . 

Let z be the empty vertex and e\ = {y±, z} , e 2 = {y 2 , z} be the edges connected to the empty 
vertex. 

The value of the generalized Parisi-Wu graph (£7), [cf. Def.(4.4)] is the same when we permute 
the edges i.e when e± = {j/2, z} and e 2 = {yi, z}. 

Acknowledgements: The authors thanks Professor Sergio Albeverio for his interest and for 
valuable discussions. 
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